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Abstract 

We study the bridge between the phenomenological mass matrix model and SO(10) 
GUT. Namely, we consider the four zero texture model in the framework of the renor- 
malizable SO(10) GUT model. This unification gives more stringent constraints than 
the case where only either model is considered. However, we can obtain good fitting 
by expanding the minimal SO(10) GUT to include 120 in addition to 10 and 126 in 
Yukawa coupling and by considering both type-I and type-II seesaw mechanisms. 
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1 Introduction 



Many papers have been published on the phenomenological mass matrix model of quark- 
lepton pp. In this approach flavour symmetry and texture zero are important ingredients. 
In most of these models, quark and lepton are discussed on the same footing. This approach 
is rather directly related with observations. However there has no base why and how quark 
and lepton are unified. The latter is the main theme of the grand unified model, which 
can say the relations between mass matrices but not on the mass matrix itself. So these 
two approaches should be complementary to each other. However, curiously enough, such 
attempts have been rather few. [2J are those of such rare attempts, in which we discussed the 
relation between four zero texture (FZT) mass matrix model and renormalizable SO(10) 
GUT [3]. Renormalizable SO(10) GUT implies that fermion-Higgs coupling is resticted in 
the renormalizable Yukawa coupling. The fundamental representation of matter multiplets 
are 16, (i=family index), and 16 x 16 = 10 + 120 + 126. So the general Yukawa coupling 
is given by 

^0,^1646,10^ + y 120iii 16 i 16 i 120 H + Y 12Gij 16 i 16 j 126 H . (1) 

In this paper we combine this model with FZT model. On the side of SO(10) we consider 
the most general model with use of both type I and type II seesaw mechanism for neutrino 
masses. 



2 Four Zero texture Quark-Lepton Mass Matrices in 
SO(10) GUT 

Phenomenological quark mass matrices have been discussed from various points of view [T]. 
In this section we review our quark and lepton mass matrix model [2]. The mass term in 
the Lagrangian is given by 

L M = -q^Muijqlj - q d Rji M dij q d LJ - l^M eij l Ljj - u' Ri M Dij u LJ 



-^L,i) c M Lij u L>j - -(v' R ,) c M RlJ v' RtJ + H.c. (2) 

with 




where M u , Ma, Md, M ej ML, and M R are the mass matrices for up quarks, down quarks, Dirac 
neutrinos, charged leptons, left-handed Majorana neutrinos, and right-handed Majorana 
neutrinos, respectively. The mass matrix of light Majorana neutrinos M v is given by 



M v = Ml — MlM R l M D , (4) 

4 In the recent papers there is a confusion on this terminology. We have called this model FZT since we 
consider up and down quark mass matrices (or lepton doublet) as a unit and they have totally 2+2=4 zeros. 
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which is constructed via the seesaw mechanism |5] from the block-diagonalization of neutrino 
mass matrix, 



M L Ml 



(5) 



M D Mr 

In the SO(10) GUT scheme, we consider the general renormalizable model, in which Yukawa 
coupling involves three types of mass matrices involving 10, 120, and 126 Higgs. These 
Higgs fields are decomposed into 

10 = (6, 1,1) + (1,2, 2), 

120 = (15, 2, 2) + (6, 3,1) + (6, 1,3) + (1,2, 2) + (10, 1,1) + (10,1,1), (6) 
126 = (10,1,3) + (10,3,1) + (15,2,2) + (6,1,1). 

under SU(A) x SU(2) L x SU{2) R . Thus 10 and 126 have two SU(2) doublets and appears 
CG coefficients 1, -3, respectively. On the other hand, 120 has four doublets and no CG 
coefficient. And 

16x16x10 D 5 H (uu c + uu c ) + 5 H (dd c + ee c ), 
16x16x120 D 5 H vv c + 45 H uu c + 5 H (dd c + ee c ) + 45 H (dd c -3ee c ), (7) 
16x16x126 D l H v c v c + 15 H vv + h H (uu c -3uv c ) + 45 H (dd c -3ee c ). 



under SU(5) decomposition. The first and second terms of 126 are the right-handed and 
left-handed Majorana masses, respectively. So six mass matrices included in SO(10) have 
the following form, 

M u = S+S'A + eS' = S U + A U , 
M d = aS + 5A + S' = S d + A d , 

M D = S + 5" A - 3eS' = S D + A D , (8) 
M e = aS + A-3S' = S e + A e , 
M L = PS' = S L , 

M R = lS' = S R . 

Here S, A, S' represent common structure for mass matrices which come from 10, 120, and 
126, respectively. a,S,5',5",e, /3, 7 are relative coefficients of vacuum expectation values 
(VEVs). Sf(Af) (/ = u,d, D,e, L, R) are the symmetric (antisymmetric) part of Mf. Fur- 
thermore,we put a ansatz that the mass matrices M u , M d , M e and M v are hermitian and 
have the same textures. That is, S a are real and A a are pure imaginary. Our model is 
different from the Fritzsch model [5] in the sense that (2,2) components are not zeros and 
that our model deals with the quark and lepton mass matrices on the same footing. The 
mass matrices Md, Ml, and Mr are, furthermore, assumed to have the same zero texture as 
M v [6]. This ansatz restricts the texture forms [2] and we choose the following our texture 
because it is most closely related with the NNI form [S]. 



NNI : * * , Our Texture :(***. (9) 
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Our model is different from the NNI model in the sense that (2,2) components are not zeros 
and that our model deals with the quark and lepton mass matrices on the same footing. 
Thus, in the FZT model, the quark and lepton mass matrices are described as follows. 



M, 



M, 



R 




M e = P e \ a e b e c e P e t = aee -ir e K ^ I (1Q) 



Mr 



where Pf = diag(e* 7/1 , e nf2 , e* 7/3 ), 77 = 7/1—7/2 and 07 = 7f2— 7/3- Hereafter we collectively 
describe these matrices as 

M f = I a e e~ iT f b f c f e ia f I (11) 






a f e lT f 





a e e~ iT f 




c f e ia f 





c f e^f 





in addition to ([9]). Here f=u,d,D,e,L,R. If we denote the mass eigen values of Mf as rrifi (i 
1, 2, 3), then a/, bf, Cf can be expressed in terms of rrifi and df. 



m fl m f2 m f3 / (df - m fl )(d f - m f2 )(df - m /3 ) 

a f - 1/ ^— • c ^ - v ^ • 

bf = mfi + mf2 + mf 3 — df. (12) 

Here 

< m/i < — m/2 < m/3 for |m/i| < df < |m/2|, 

< — m/i < m/2 < m/ 3 for |m/ 2 | < df < |m/ 3 |. (13) 

The number of parameters of quark and charged lepton are from f=u,d,e of (151) , 

S,A,S' ^4 + 2 + 4= 10 (14) 

and 4 coefficients of VEVs, a, 5, 5', e, and totally 10 + 4 = 14. On the other hand, the number 
of constraints from experiments are quark masses (6), CKM mixing angles (3), the Dirac 
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phase (1) and the lepton masses (3), that is, 6 + 3 + 1 + 3 = 13. Therefore only 1 parameter 
is remained after the data fitting of quark and charged lepton. Using this one free parameter 
and 5", j3, 7, we can fix light neutrino masses (3) and MNS mixing angles (3) and the phases 
(3). Among these the direct experimental constraints are A 2 m 23 ,A 2 m 12 and three MNS 
mixing angles. 



3 Numerical analyses 

FZT matrices are diagonalized as 

UjM f U f = diag(m/i, rap, m /3 ) 

Here 



U f = PjOf, P f = diag(l, 77, o s + T f ) = (1, a/2, a 



and 



O, 



where 



I I (df-mf 1 )mf 2 mfs 

R fl d f 



(d f -m fl )m fl 



Rfi 



(df-mf 2 )rnf 1 mf 3 

(df-mj2)mf 2 
Rf2 



(df-m.f A )mf 2 mfi \ 

R fi d f 
(d f ~m f3 )m f3 



m fl (d f -m f2 )(d f -m f3 ) 
R fl d f 



mf 2 (d f -m f3 )(d f -m fl ) / m f3 (d f -m fl )(d f -m f2 ) 



(15) 
(16) 



(17) 



R f i = (m fl - m f2 ){m fl - m /3 ), R f2 = (m /2 - m /3 )(m /2 - m fl ), 

Rf3 = {m f3 -m fl ){m f3 -m f2 ). (18) 

Hereafter we denote (m/i, m/ 2 , m/ 3 ) for f = u,d, and e as (m u ,m c ,m t ), (m d ,m s ,mb), and 
(m e , m,,, m r ), respectively. The CKM quark mixing matrix Uqkm = U}jJ d is given by 



CKM ) 12 



- l^x u x d - e^/^(l - Xu )(l- Xd ) 



m. 



\m d \m s 



CKM) 23 



(Uc Km) 13 ~ 

and the Dirac phase is 



e ** W ™x u (l - x d ) + e^W^(l _ Xu)Xd . 

/ J U — 2^ + e iQ2 Vx u (l - x d - e lQ3 V I - x u ). 

y m c m l hXd 



(19) 



5 9 ~ arj 



(e iQ V(l -a; u )(l-x d ) +e iQ2 v ^^)* 

(e iQ 3 ^(1 - ijirf - e™* ^(1 - x d )) (e^ y/(l- Xu )x d - e ia * y/x u {l-x d ))* ' 



(20) 
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for \{Uckm)\z\ *C 1, where we define 

df 

x f = 

m f3 

«2 = a U 2 - 0L d2 = r u -T d = At, (21) 
«3 = ««3 - a d3 = r u -r d + (a u - a d ) = At + Act. 

We note that the values of m;, affects sensitively on the CKM quark mixings, but the value 
of mt does not. 

There are additional constraints among the parameters in the FZT model, when it is 
embedded in the SO (10) GUT model. Let us first discuss the consraints in the quarks and 
charged leptons sector. In the SO(10) frame mentioned above, Sf and Af (/ = u, d, e) satisfy 
following relations: 

4aS u = (3 + ae)S d +(l-ae)S e , (22) 
SA U = 6'A d = 66'A e . (23) 

Expressing in terms of the components, we have the following constraints among the param- 
eters: 

Aaa u cos (At + T d ) = (3 + ae)a d cos T d + (1 — ae)a e cos r e , (24) 

4ac u cos(Acx + a d ) = (3 + ae)c d cos a d + (1 — ae)c e cos cr e , (25) 

Aab u = (3 + ae)b d + (1 - ae)b e , (26) 

Aad u = (3 + ae)d d + (1 — ae)d e , (27) 

Sa u sin( Ar + T d ) = 5'a d sin T d = 5S'a e sin r e , (28) 

Sc u sin(Acr + a d ) = 5'c d sin a d = SS'c e sin a e . (29) 

From these constraints we have 

F(r, d u , d d ) 2 [Aaa u cos(Ar + T d ) - (2 + K,)a d cosT d } 2 - [4ac u cos(Aa + a d ) - (3 + n)c d cos a d ] 2 

= (l-K) 2 [a 2 e F(r 1 d UJ d d ) 2 -c 2 e ] 1 (30) 



with constraints 



4aa u cos(Ar + r d ) - (3 + n)a d cosT d 

-l<COST e = < 1, (31) 

(1 - K)a e 

4ac u cos(Aer + a d ) - (3 + K)c d cosa d 

-1 < coscr e = ; s < 1- 32 

(l-K)Ce " 



where 



^ „/ , , v Qsma d c e sin(7 e c e sma e 

n = ae, r = — , F(r, d u , of d ) = : = : = : 6. (33) 

6 a d smr d a e smr e a d smr d 
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We note that r^ and ad are obtained from Eqs. fl24|) - fl29|) as 

a u sin Ar c u sin Act 

tanr d = — , tana d = — . (34) 

rcid — a u cos Ar red — c u cos Aa 

And a and e are determined as functions of d u , dd, and d e : 

a (m d + m s + m b )d e - (m e + m M + m T )d d 

(m d + m s + m h - m e - m M - m r )c/ u - (m u + m c + m t )(d d - d e ) ' 
^ _ K + ra c + m t )(3d d + 4) - [3(m rf + m A + m 6 ) + (m e + m M + m T )] d M 
(m d + m s + m b - m e - - m T )d u - (m u + m c + m t )(d d - d e ) 



Note also that from Eq. f[33|) . when | cosr e | = 1, | coscr e | becomes 1, Therefore, in fact, we 
need only one of Eqs. ( 13~TT) and ( |32l ) . 

As was shown in Eqs. ( 1T91) and ( 120|) . the parameters d u , dd, Ar, and Aa are constrained 
by the observed CKM mixing matrix. The best fit is realized for the following values of the 
parameters, 

Ar = tt/2, (37) 
Aa = -0.121, (38) 
d u = 0.9560 m t , (39) 
d d = 0.9477 m b , (40) 

which are used in the following analysis. Then r^ and ad are determined by two parameter 
d e and r from Eqs. (134"|) since the parameters aj, bf, Cf are fixed by df from Eq. (fl~2l) . For 
these best-fit-parameters mentioned above, we obtain 

(41) 
(42) 
(43) 
(44) 

Here we have used the best fit values of the following quark and charged lepton masses 
estimated [9] at the unification scale /i = Mx- 

1.33i°i£MeV, 
26.51^ MeV, 

1.00 ± 0.04 GeV, (45) 
68.59813 ± 0.00022 GeV, 



Therefore we have two parameters d e and r left to be fixed, which are determined from 
Eq. fl30|) with Eqs. fl3Tl) and fl32l) . The solution of these equations is depicted in Figure 1 



\{UcKM, 


1 12 1 


= 0.2251, 


\{UcKM, 


*23 1 


= 0.0340, 


\{UcKM, 


1 13 1 


= 0.0032, 




6 q 


= 58.86°. 



m u (M x ) 


= 1.041^ MeV, 


m d {M x 


m c (M x ) 


= 302^5 MeV, 


m s (M x 


m t {M x ) 


= 129+^ 6 GeV, 


m b (M x 


m e {M x ) 


= 0.32502032 ± 0.00000009 MeV, 


m M (M x 


m T {M x ) 


= 1171.4 ±0.2 MeV. 
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in the d e — r plane. As seen in Fig. 1 there are two allowed curves which we call Sol. (a) 
and Sol. (b), respectively. The typical values of d e /m T for these solutions are 0.935883 for 
Sol. (a) and 0.307197 for Sol. (b), which we shall use in the following analysis. 

Thus we have succeeded to fit 13 parameters consistently in quark and charged lepton 
sectors. Using these parameters and M D , M R and M L , now let us proceed to discuss the 
neutrino masses and lepton mixings. As seen from Eq. (Q we have three more free param- 
eters 6", (3, and 7 in the neutrino sector. The neutrino mass matrix is given by the seesaw 
mechanism as 

M y = M L - M D M^Ml. (46) 

The diagonalization of the charged lepton mass matrix M e is simmilar to the quark mass 
matrix. On the other hand, since the mass matrix for the Majorana neutrinos is symmetric, 
M v is diagonalized as 

UlMJJt = diag (|m x I, m 2 ,m 3 ) , (47) 

where |m 1 |,m 2 ,and m 3 are real positive neutrino masses and the unitary matrix U v is de- 
scribed as 

U v = P}O v Q v . (48) 

Here, in order to make the neutrino masses for the first generation to be real positive, we 
introduce an additional diagonal phase matrix Q v defined by 

Q^diagfrM). (49) 

The MNS lepton mixing matrix Umns of the model is given by 

Umns = U\JJ V = OjP £ O u Q u , (50) 

where P^ = P e Pl is diagonal phase matrix. Eq. (!50|) is changed to the standard representation 
of the MNS lepton mixing matrix as well as the CKM quark mixing matrix, 

f~X fX f~X rX rX 

c 13 c 12 c 13 i 12 "13 

s~X q£ q^> s~X' qI fD^^X (~X /-X q^> £3^^\ S~X' 

c 23' b 12 i, 23 L 12' b 13 e c 23 c 12 " b 23 5 12 5 13 e 6 23 c 13 

Q I Q I f~X /-X fD^^l O ^ S~X S~X' Q I SD^^l S~X (~X 

*23*12 L 23 c 12' b 13 e 5 23 c 12 °23 6 12' b 13 e C 23 L 13 

xdiag(l,e^ 2 ,e i03 ). (51) 

Here Q comes from the rephasing in the charged-lepton fields, 5i is the Dirac phase, and fa 
i — 1,2 are the Majorana phases in the MNS lepton mixing matrix. 

From Eq. (157J1) the neutrino oscillation angles and phases of the model are related to the 
elements as follows: 

t 2 a \(Umns)i2\ 2 fKO , 

lan (7 so i ar — j-— r rr- , Ky^i 



|2 ' 

sin a 20 atm = 4\(Umns)23\ 2 \(U MN s)33\ 2 ■ (53) 



\{Umi\ 

,2 , 
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In Fig. 2 - Fig. 5, we present the numerical values of the mixing angles in terms of the 
parameter f3j. Here we have used the following experimental constraints from global analysis 
of neutrino parameters |10j . 

0.25 < sin 2 fl 12 < 0.38, (54) 
0.35 < sin 2 # 23 < 0.65, (55) 
sin 2 # 13 < 0.03. (56) 

and 

Am 2 21 = (7.2-8.9) x 10~ 5 [eU 2 ], (57) 
Am 2 2 = (2.1-3.1) x 10" 3 [e1/ 2 ]. (58) 

Then the neutrino mass square ratio becomes 

Am 2 21 /Am 2 32 = 0.023 - 0.042. (59) 

As is seen in Fig. 2 - Fig. 5 (especially Fig. 4) we have consistent values of the lepton mixing 
angles with the observed experimental data if we take 5" = and tune (3^. Namely we 
obtain 

sin 2 2 = 0.32, sin 2 3 = 0.52, sin 2 3 = 3.0 x 10" 6 , for f3j ~ 3.4515 x 10 5 at Sol. (a), (60) 
sin 2 2 = 0.32, sin 2 3 = 0.53, sin 2 3 = 2.3 x 10" 4 , for /3 7 ~ 2.9240 x 10 7 at Sol. (b). (61) 

On the other hand, for these parameters, the model predicts somehow small value for ^ m jf , 

Z\m 32 

as seen in Fig. 6. 

— -|i = 1.7 x 10~ 7 for Sol. (a), — -f- = 3.4 x 10" 4 for Sol. (b). (62) 



4 Discussions 



In this paper we have considered the bridge between the phenomenological FZT model and 
SO(10) GUT. Namely we incorporated FZT in the framework of SO(10) GUT. So our model 
suffers more stringent constraints than tha case when either only FZT or only SO(10) GUT 
was considered. Neverthless, our model gave values consistent with the observations of 
quark-lepton. Only exception is the mass square ratio of neutrino. However, this is not 
so bad news. For the other observables' fittings are very nice, and FZT model is itself an 
approximation. Also we have not considered the renormalization group equation (RGE) for 
neutrino sector. Namely we bottomed up the low energy spectrum to GUT and fixed the 
Mfl and could construct light neutrino mass matrix first at GUT. So we must top down this 
matrix to the electroweak scale and fit with oscillation data [11] [12] at this scale. However, 
RGE effect works little on neutrino mixing angles and mass ratios if neutrino masses are not 
degenerate. Rather they depend on RGE of quark masses in which several ambiguities are 
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m t x 2 


mt/2 


m b x 2 


m b /2 


m T x 2 


m T /2 


Pi 


1.15 x 10 8 


7.30 x 10 6 


8.61 x 10 7 


7.60 x 10 4 


1.24 x 10 6 


2.95 x 10 4 


sin 2 #12 


0.315 


0.333 


0.311 


0.330 


0.312 


0.315 


sin 2 # 23 


0.573 


0.481 


0.582 


0.460 


0.917 


0.131 


sin 2 9 13 


2.15 x 10" 4 


3.02 x 10" 4 


2.30 x 10" 4 


2.20 x 10" 4 


1.55 x 10~ 3 


2.25 x 10" 4 


Am 2 2 /Am 2 3 


8.79 x 10" 4 


2.87 x 10" 4 


9.57 x 10" 4 


2.60 x 10" 4 


1.07 x 10" 7 


1.77 x 10" 4 



Table 1: The values around the peak when changing the third generaton masses. 



left unsolved. Nobody knows clearly what is happening above TeV scale. In most cases, the 
ratio mfi/mf2 is stable but m/i/ra/ 3 is not. Therefore, by changing only m/3 in the very 
wide range we will be able to check the stability of the bi-large mixing. 

We note that we must check carefully when we change mj because the CKM quark mixing 
matrix is sensitive to mj. Therefore only when we change the mass of bottom quark, we use 
the following values: 

• In the case m b = (l.GeV) x 1/2, 

At = tt/2, Aa = -0.314, d u = 0.993mf, d d = 0.993m 6 , (63) 



\{U C km)i2\ = 0.2254, |([W) 23 | = 0.0347, \{U CKM )iz\ = 0.0031, 5 q = 60.67°. 

(64) 

• In the case m b = (l.GeV) x 2, 

At = tt/2, Aa = -0.063, d u = 0.811m t , d d = 0.796m b , (65) 



\{U C km)i2\ = 0.2245, \(U C kmU\ = 0.0327, |(fWr) 13 | = 0.0032, 5 q = 59.32°. 

(66) 

We show, in Table 1 and Fig. 7, how the lepton mixing angles and neutrino masses are 
sensitive to the quark and charged lepton masses of the third generation by changing factor 
2 of them. We find that RGE effects to our solution of the bi-large lepton mixing are small 
from m t and m T , while large from m&. 
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Figure 1: The black lines satisfies Eq. (1301) . The blue and red regions show the additional 
conditions Eqs. (l3Tj) and (|32|) . respectively. The two overlapping regions of black line, red 
and blue areas in (i) are the allowed regions. The panels (ii) and (iii) are the enlarged figures 
of the allowed regions of (i) and are called Sol. (a) and Sol. (b), respectively, hereafter. 
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Figure 2: The Solution (a) in normal hierarchy. The red, blue, violet areas satisfy the solar 
Eq. (60), atmospheric Eq. (61), and CHOOZ Eq. (62), respectively. The green area shows 
the mass square ratio Eq. (65). (i) The diagram of /?7 versus 5" . (ii) The enlarged diagram 
of (i). (iii) Mixing angles as functions of ^7: The solid red, blue, and violet lines depict 
sin 2 6'i2, sin 2 6> 2 3, and sin 2 ^, respectively. 



Figure 3: Same as Figure 2 but Sol. (a) in inverted hierarchy. There is no red region in (i) 
and (ii) because sin 2 # 12 ~ 1. 



Figure 4: Same as Figure 2 but Sol. (b) in normal hierarchy. There is no violet region in (i) 
and (ii) because sin 2 # 13 <C 1. 



Figure 5: Same as Figure 2 but Sol. (b) in inverted hierarchy. There is no red region in (i) 
and (ii) because sin 2 #12 — 1. 



Figure 6: The neutrino masses in Sol. (a) and (b). The is the mass eigenvalues of light 
neutrinos. The red, blue, and green lines represent the neutrino masses in light order. 



Figure 7: The third generation mass dependence of the mixing angles as function of ^7. 
The red, blue and violet lines show sin 2 6' 12 , sin 2 #23 and sin 2 # 13 , respectively in Sol. (b) in 
normal hierarcy. The top panel shows the values when we use the best values of masses in 
Eq. (145]) . Other panels show whether the bi-large mixing angle is reproduced or not if the 
third generation mass would take double or half value. Sol. (a) has an behaviour similar to 
Sol. (b) 
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